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FOREWORD 


This  is  the  Summary  Volume  for  three  separately  bound  volumes  reporting  research 
completed  under  the  general  terns  of  the  Office  of  Civil  Defense  Subtask  Number  4113E, 
"Sensitivity  Analysis  of  Civil  Defense  Systems  and  Components."  Volume  I  by  Floyd 
Guess  reports  on  A  Cost/Effectiveness  Computer  Procedure  for  Optimum  Allocation  of 
Fallout  Shelter  System  Funds  Under  Uniform  or  Variable  Risk  Assumptions.  Volume  II 
by  Rodney  Sink  reports  on  A  Sensitivity  Analysis  of  Selected  Parameters  Based  on 
8  SHSA's.  Volume  III  by  John  Neblett  and  K.E.  Willis  reports  on  A  Generalized 
Sensitivity  Analysis  of  CD  Systems. 

Hie  research  of  che  author*}  was  very  ably  supported  by  Herbert  Hill,  Helen 
Anderson,  and  Mary  Woodside.  Philip  McMullan  and  Robert  Brooks  are  acknowledged 
for  their  valuable  assistance  in  the  preparation  of  the  final  report. 


iii 


ABSTRACT 


This  document  summarizes  a  three  part  study  concerning  sensitivity  analysis  of 
CD  systems  and  components  in  a  fallout  environment.  In  the  first,  a  cost/eff ective- 
ness  computer  program  is  developed  for  optimum  allocation  of  fallout  shelter  system 
development  funds  under  uniform  or  variable  risk  assumptions.  This  program,  intended 
for  use  in  OCD  planning  studies,  is  programmed  for  the  CDC  3600.  It  is  applied  in 
example  studies  using  data  on  OCD  Region  6.  The  second  part  of  the  study  is  a  sen¬ 
sitivity  analysis  of  selected  parameters  based  on  8  SMSA's.  It  employs  the  transpor¬ 
tation  algorithm  in  a  study  of  movement  of  people  to  fallout  shelters.  The  results 
show  how  estimated  casualties  vary  as  movement-to-shel ter  patterns  vary  from  restric¬ 
tion  to  a  standard  location  up  to  frte  movement  within  the  SMSA.  They  also  Indicate 
that  detailed  planning  for  shelter  utilization  can  be  very  effective  in  reducing 
expected  fallout  casualties  when  the  number  of  shelter  spaces  exceeds  the  population 
of  an  SMSA.  In  the  third  part  of  the  study,  a  generalized  sensitivity  analysis  is 
made  of  the  parameters  used  in  fallout  vulnerability  analysis  models  which  determine 
total  dose  and  equivalent  residual  dose.  The  results  show  that  fallout  casualties 
are  very  sensitive  to  estimates  of  initial  intensity,  protection  factor,  and  fallout 
decay  rate;  and  the  equivalent  residual  dose  during  stay  in  the  initial  shelter  is 
insensitive  to  the  recovery  fraction  and  rate  of  recovery  parameters  of  the  equivalent 
residual  dose  equation. 
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Introduction  and  Summary  to  Sensitivity  Analysis 
of  Civil  Defense  Systems  and  Com>onent» 


I.  INTRODUCTION 


A.  Shelter  Systems  Analysis 

The  general  Bubject  of  shelter  systems  analysis  is  outlined  in  this  section  in 
order  to  place  in  context  the  sensitivity  analysis  reported  in  the  following 
sections . 

The  existing  Civil  Defense  system  consists  principally  of  fallout  shelters, 
their  support  «nd  operational  plans  for  tbclr  use.  Information  from  the 

National  Fallout  Shelter  Survey  (NFSS)  gives  the  numbers  and  locations  of  spaces 
that  could  be  used  and  their  respective  protection  factors.  Information  also  exists 
as  to  supplies  and  equipment  contained  in  these  shelters.  Operational  plans  of 
varying  degrees  of  completeness  also  exist  throughout  the  country  for  the  use  of 
shelters  in  an  emergency- -that  is,  for  population  assignment  to  these  spaces  and 
shelter  management  policies. 

The  evaluation,  or  vulnerability  analysis,  of  the  shelter  system  requires  esti¬ 
mation  of  both  the  ef fectiveness  of  the  system  as  it  now  exists  in  alleviating  the 
effects  of  nuclear  attack  and  the  effectiveness  of  alternative  systems. 

In  order  to  perform  an  analysis  of  the  effectiveness  of  an  existing  or  proposed 
shelter  system,  it  is  necessary  to  consider  three  problem  areas.  First,  some  assump¬ 
tions  must  be  made  about  the  radlation--its  initial  intensity,  the  time  it  will  ar¬ 
rive,  the  time  span  during  which  it  will  continue  to  fall,  and  the  rate  of  decay  o£ 
radioactivity.  Second,  the  shelter  posture  must  be  defined,  principally  In  terms 
of  the  protective  characteristics  of  the  shelter  and  the  time  required  to  reach 
shelter.  Third,  some  assumptions  must  be  made  about  the  physiological  effects  re¬ 
sulting  from  the  assumed  levels  of  radioactivity. 

Systems  analysis  also  requires  that  measures  of  effectiveness  be  established 
so  that  the  effect  of  alternatives  can  be  compared.  Some  coianon  measures  of  effec¬ 
tiveness  are  probability  of  fatality,  probability  of  non-fatal  casualty,  or  re¬ 
quired  stay  in  shelter  to  meet  some  specified  criterion--such  as  equivalent  residual 


dose  (ERD)  returning  to  an  allowable  level. 

The  results  of  systems  analysis  are  used  to  evaluate  the  current  civil  defense 
system,  co  evaluate  alternative  systems,  and  hence,  to  insure  that  the  limited  funds 
for  system  improvement  yield  the  greatest  gains  possible. 

Assumptions  about  thermonuclear  war  environments  obviously  are  questionable  be¬ 
cause  of  the  lack  of  observational  experience.  Uncertainty  also  results  from  lack 
of  reliable  data  for  some  necessary  inputs  to  the  analysis,  such  as  the  dose-response 
curve  for  fallout  casualties.  To  evaluate  the  magnitude  of  the  uncertainty  in  re¬ 
sults  and  conclusions,  one  must  test  the  effects  of  uncertainties  in  these  input 
variables.  W^«n  these  effects  are  evaluated,  the  degree  of  confidence  which  can  be 
placed  in  the  results  of  systems  analysis  will  be  better  known;  further,  research 
can  be  directed  in  some  cases  towards  reducing  the  uncertainty  of  the  inputs  to 
which  evaluations  are  most  sensitive.  A  collateral  benefit  for  systems  analysis  is 
simplification  of  the  treatment  of  variables  to  which  the  outputs  are  insensitive. 

In  some  cases,  the  less  sensitive  variables  may  be  eliminated. 

The  sensitivity  analysis  carried  out  during  the  past  year  is  described  next. 


B.  Objectives  of  the  Study 

The  contract  language  is  quoted  below: 

"Using  the  tools  developed  under  OCD  Projects  4631A  and 
4521A  (formerly  4104A  and  4104B  respectively),  perform  for 
selected  local  areas  detailed  sensitivity  analyses  of  the  variable 
elements  of  shelter  systems  evaluation  (such  as  relative  location 
of  people  and  shelter,  warnings  and  reaction  time,  movement,  to 
shelter,  shelter  stay  time,  food  and  water  supplies,  mass  decontami¬ 
nation  and  evacuation)  over  the  range  of  feasible  fallout  environ¬ 
ments.  Emphasis  would  be  placed  on  incorporating  the  cost-effec¬ 
tiveness  findings  of  concurrent  research  on  shelter  and  support 
systems . 

"Use  the  results  of  such  analysis  to  (1)  assist  in  the 
design  and  performance  of  total  CD  system  and  (2)  assist,  as 
specifically  authorized  by  separate  letter  instruction,  in  the 
formulation  of  computer  programs  and  systems  evaluations  using 
the  NREC  and  OCD  damage  assessment  system  for  the  1103A  or  CDC 
3600  computer." 

The  description  of  the  approach  and  research  findings  is  presented  in  three 
volumes  of  this  report,  as  follows: 
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Volume  No. 


Title 


I  A  Cost/Effectiveness  Computer  Procedure 
for  Optimum  Allocation  of  Fallout  Shelter 
System  Funds  Under  Uniform  or  Variable 
Risk  Assumptions 

II  A  Sensitivity  Analysis  of  Selected  Parameters 
Based  on  8  SMSA's 

III  A  Generalized  Sensitivity  Analysis  of  CD 
Systems 

The  summary  below  follows  the  same  sequence.  Sensitivity  analysis  of  various 
components  of  CD  systems  is  the  subject  of  Volumes  II  and  III,  while  cost  and 
effectiveness  analysis  is  the  subject  of  Volume  I.  The  marginal  expenditure  pro¬ 
cedure  described  in  Volume  I  provides  for  optimal  allocation  of  funds  to  shelter 
improvement  options,  subject  as  always  to  certain  policy  decisions. 

All  elements  of  the  research  were  directed  to  the  “fallout  only"  attack  environ 
ment.  Reconsnendations  are  offered  for  future  research  to  incorporate  prompt  weapon 
effects;  a  small  effort  also  will  be  necessary  to  improve  the  models  and  methods 
employed . 

The  next  three  sections  summarize  the  research  approach,  relevant  findings,  con 
elusions,  and  recommendations  for  each  of  the  three  principal  tasks  and  chapters. 

II.  SUMMARY 

A.  A  Cost/Effectiveness  Computer  Procedure  for  Optimum  Allocation  of  Fallout 

Shelter  System  Funds  Under  Uniform  or  Variable  Risk  Assumptions 

1.  Approach  to  the  Research 

The  dynamics  of  civil  defense  planning  and  systems  evaluation  require  a 

procedure  that  will  readily  yield  approximate  answers  to  such  questions  ac : 

a.  How  many  lives  could  be  saved  with  an  expenditure  o.e  $1  billion, 

$5  billion,  $15  billion,  or  any  other  arbitrary  budget  level? 

b.  How  should  such  budgets  be  exp ended --with in  regions,  states, 
metropolitan  areas,  counties,  and  smaller  locations — for  each 
assumed  budget  level? 

c.  How  and  where  would  the  expenditure  pattern  vary  if  one  assumed  that 
all  parts  of  the  United  States  were  subject  to  a  uniform  hazard 
level--or  to  variable  hazard  levels? 
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d.  What  la  the  dollar-coatlng  relationship  between  attack  hazard, 
casualties,  and  standards  of  acceptable  shelter? 

Our  approach  has  been  to  develop  and  demonstrate  a  computerized  procedure 
that  will  answer  these  and  related  questions  under  the  following  constraints: 

a.  Use  available  data,  such  as  the  NFSS  and  National  Location  Code, 
wherever  feasible. 

b.  Use  existing  computers,  programs,  and  planning  techniques  as  used 
by  OCD,  such  as  RISK  or  NAHICUS ,  wherever  feasible. 

c.  Provide  means  for  readily  accepting  new  or  improved  data,  and  for 
using  interim  estimates  where  critical  data  deficiencies  exist. 

d.  Concentrate  initially  on  cost  and  effectiveness  of  a  fallout  shelter 
system,  but  Include  provision  for  extending  the  system  to  blast  and 
other  prompt  effects  considerations. 

e.  Provide  means  for  acceptipg  data  on  all  civil  defense  subsystems  of 
measurable  significance. 

At  this  point  in  the  research,  it  is  prudent  to  say  only  that  we  present 
a  rational,  well-defined  methodology  for  engaging  in  '’risk-oriented  program¬ 
ming";  it  is  not  part  of  this  research  task  to  examine  the  advisability  of 
risk-oriented  courses  of  action.  Should  "risk-oriented  programming"  be  decided 
upon,  the  tool  developed  in  Volume  I  is  ready  for  application.  Further  analysis 
can  offer  guidance  in  making  the  decision  of  whether  or  not  to  adopt  "risk". 

In  order  to  apply  the  risk-oriented  budget  allocation  procedure,  the 
following  input  data  are  required  (note  that  this  demonstration  is  limited  to 
the  fallout-only  case)  : 

a.  The  present  shelter  posture  and  methods  for  estimating  casualties 

(or  other  measure  of  effectiveness)  must  be  specified.  This  requires 
a  population  data  base  (here  taken  as  residential),  a  means  of 
assigning  them  to  shelter  (of  a  specified  protection  factor),  and  a 
casualty  criterion. 

b.  An  estimate  is  required  of  the  probable  attack  environment  defined 
by  areas  of  such  a  size  that  they  experience  essentially  a  homogene¬ 
ous  environment.  For  fallout  only,  such  areas  may  be  quite  large. 

A  further  constraint  on  the  size  of  the  area  within  which  one  matches 
shelter  and  population  is  that  the  two  are  close  enough  together  that 
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shelters  may  be  assumed  occupied  on  arrival  of  fallout.  The  accumu¬ 
lation  unit  used  here  for  shelter  and  population  is  the  county.  For 
illustrating  the  method,  NREC  RISK-type  attack  environments  are 
assumed  at  two  levels:  50%  (expecced  value)  and  95%  (only  1  in  20 
chances  of  being  exceeded). 

c.  Further,  a  specified  set  of  alternatives  for  spending  for  shelter 
improvement  and  a  means  for  computing  the  contribution  of  each 
alternative  to  the  measure  of  effectiveness  are  required.  The 
alternatives  chosen  in  this  analysis  are  based  on  NFSS  Phase  2  esti¬ 
mates  of  improvable  shelter  (improvable  by  shielding  or  ventilation) 
and  PF  Category  1  shelter  from  NFSS  Phase  1  data.  When  existing 
and  improvable  NFSS  shelters  are  exhausted,  new  shelter  spaces — for 
not  over  10%  of  the  population-- are  considered  to  be  available  at 
$25  per  space.  For  costing  complete  shelter  coverage,  an  additional 
alternative  of  unlimited  new  construction  at  $50  per  space  is  also 
examined.  With  these  input  data,  the  appropriate  measure  of  effec¬ 
tiveness  (survivors  t.dded,  for  example)  may  be  calculated  for  combi¬ 
nations  of  spending  alternatives.  If  one  chooses  a  measure  of  effec¬ 
tiveness  having  a  continuum  of  values  (such  as  probability  of 
fatality,  from  which  survivors  added  are  derived),  the  optimum  fund 
allocation  problem  is  one  of  classical  linear  programming.  This 

may  be  solved  by  well-known  methods,  but  a  program  of  realistic  size 
can  be  very  time  consuming  even  when  a  large-scale  computer  is  used. 

On  the  other  hand,  as  done  here,  if  one  chooses  an  objective  function 
which  takes  on  two  discrete  values,  such  as  "ERD  greater  than  200r 
is  ineffective",  the  optimal  allocation  problem  reduces  to  one  of 
ordering  the  shelter  improvement  alternatives  in  the  decreasing  order 
of  cost/effectiveness.  Then  one  simply  chooses  the  shelter  improve¬ 
ment  alternative  of  highest  cost/effectiveness ,  then  the  one  with 
the  next  highest,  and  so  on,  until  either  the  budget  or  the  spending 
alternatives  are  exhausted. 

d.  A  finite  or  infinite  budget  must  be  specified.  Depending  upon  the 
budget  and  the  particular  improvement  alternatives  specified,  one 
may  "buy"  all  the  available  alternatives;  more  commonly  in  practice, 
the  budget  is  the  constraint  and  not  all  alternatives  will  be  "bought". 
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2.  Demonstration  of  the  Model 

To  demonstrate  the  model,  an  analysis  is  made  of  OCD  Region  6.  The  total 
population  of  these  8  states  (619  counties)  is  14,000,000.  In  all,  14  cost 
studies  are  run,  using  selected  combinations  of  the  following  parameter  values: 

a.  The  budget  level  (unlimited,  $20  million,  and  $10  million); 

b.  The  risk  level  (50%,  or  "expected  value"  of  attack  environment;  and 
95%,  a  very  "safe"  planning  criterion); 

c.  The  Maximum  ERD  defining  an  acceptable  shelter  Protection  Factor- 
fallout  level  combination  (lOOr,  threshold  of  perceptible  physiological 
effects;  and  200r,  threshold  of  casualties); 

d.  Residential  basements  excluded,  or  included; 

e.  Residential  basements  and  Category  1  NFSS  shelters  excluded; 

f.  Shelter  overcrowding  by  25%; 

g.  Ventilation  costs  reduced  to  $3  per  space;  and 

h.  PF's  of  shelter  multiplied  by  approximately  2  (reflecting  approxima¬ 
tions  of  the  order  of  magnitude  of  the  conservative  bias  used  in  the 
Phase  1  NFSS  calculations  when  compared  to  the  more  nearly  correct 
calculations  using  the  current  Engineering  Manual  data). 

The  program  is  designed  to  make  provision  for  prompt  weapons  effects,  but 
they  are  temporarily  excluded  from  the  analysis,  Their  later  inclusion  is 
discussed  in  Volume  I. 

3.  Conclus ions 

a.  The  demonstration  applied  to  Region  6  shows  the  practicality  of  carrying 
out  large-scale  marginal  cost/effectiveness  analyses. 

b.  There  is  a  great  need  for  reliable  input  data--particularly  the  unit 
cost  and  available  numbers  of  new  shelters. 

c.  The  model  and  case  studies  described  give  much  guidance  about  shelter 
development  planning.  However,  the  existing  procedure  still  leaves  scope 
for  the  decision  maker.  His  willingness  to  apply  this  procedure  is  con¬ 
tingent  on  the  sensitivity  of  allocation  decisions  to  judgment  factors  such 
as  listed  below.  Model  case  studies  can  give  valuable  insights  concerning 
these  questions. 
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(1)  Should  "risk-oriented"  programs  be  developed? 

(2)  What  will  be  the  severity  of  the  impact  of  future  changes  in 
expected  attack  environments? 

(3)  What  compromise  plan  would  be  acceptable  for  an  area  in  which 
the  optimum  improvement  programs  are  very  different  for  day 
and  night  population? 

4.  Recommendations 

a.  It  is  recommended  that  the  marginal  cost/effectiveness  model  be  ex¬ 
tended  to  incorporate  prompt  weapons  effects,  and  that  alternative 
measures  of  effectiveness  be  explored. 

b.  It  is  recommended  that  numerous  systematic  case  studies  be  performed 
to  test  the  impact  of  possible  wrong  decisions  (due  primarily  to  the  in¬ 
herent  uncertainty  of  forecasting  the  future). 

B.  A  Sensitivity  Analysis  of  Selected  Parameters  Based  on  8  SMSA's 
1.  Approach  to  the  Research 

The  overall  objective  of  this  phase  is  the  assessment  of  casualties  in 
8  representative  cities  over  a  range  of  fallout  environments  and  shelter  utili¬ 
zation  patterns.  The  cities  are  selected  on  the  basis  of  "judgment  sampling" 
after  making  a  two-way  classification  according  to  population  density  (residents 
per  square  mile)  and  shelter  coverage  of  the  residential  population.  The 
population  density  ranges  from  1320  for  Mobile,  Alabama  to  12,950  for  Charleston, 
South  Carolina.  Shelter  coverage  (Category  1  or  better,  based  on  Phase  1  NFSS) 
ranges  from  207.  for  Orlando,  Florida  to  1607.  for  Richmond,  Virginia.  Residential 
population  ranges  from  74,000  for  Pittsfield,  Massachusetts  to  408,000  for 
Richmond,  Virginia.  The  remaining  three  cities  are  Tucson,  Arizona;  Macon, 
Georgia;  and  Waterloo,  Iowa.  The  fallout  environments  used  range  from  a 
reference  intensity  of  600  r/hr  with  time  of  arrival  of  7  hours  to  a  reference 
intensity  of  30,000  r/hr  with  time  of  arrival  of  1  hour.  The  former  approximates 
the  median  environment  based  on  an  assessment  of  prob<r  ',e  fallout  environments; 
e.g.,  about  half  are  more  intense.  The  latter  reference  intensity  is  rarely 
exceeded  (less  than  one  percent  are  more  intense). 

Shelter  assignments  are  made  using  several  methods,  including  the  well- 
known  transportation  algorithm  which  minimizes  the  expected  casualties  while 
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allowing  for  movement  in  fallout.  The  cost  matrix  for  the  transportation 
algorithm,  when  computed  (using  an  RTI  computer  program),  shows  the  Maximum 
ERD  acquired  for  each  pairing  of  population  in  a  standard  location  (SL)  and  each 
shelter  (taking  account  for  its  protection  factor).  Distances  traveled  are 
rectangular  distances  between  centers  of  SL's.  Speed  is  always  a  constant, 
conservative  2  mph.  The  fraction  of  casualties  then  is  computed  using  a  "go/ 
no-go"  definition  of  a  casualty  (ERD  >  200r). 

Since  a  shelter  assignment  plan  based  on  the  transportation  algorithm  is 
mathematically  the  best  possible  for  a  given  fallout  environment.  And  instanta¬ 
neously  takes  account  of  the  particular  environment  encountered,  it  may  be  re¬ 
garded  as  unrealistic.  (This  is  not  to  say  that  a  more  workable  decision  rule 
for  shelter  assignment  would  not  lead  to  a  similar  assignment  and  equivalent 
movement  pattern.)  Accordingly,  for  all  cities,  alternative  shelter  assignments 
are  made,  introducing  more  realistic  movement  restrictions; 

a.  Movement  restricted  to  own  standard  location. 

b.  Movement  restricted  to  within  2  miles  of  own  standard  location. 

c.  Movement  unrestricted  within  the  Standard  Metropolitan  Statistical 

Area  (SMSA). 

In  each  of  these  three  cases,  shelters  are  assigned  using  a  "hand"  map 
technique.  Dose  due  to  movement  in  fallout  is  not  computed,  as  the  circum¬ 
stances  are  assumed  to  be  equivalent  to  a  long  tactical  warning  or  to  strategic 
warning  which  would  allow  movement  to  shelter  before  arrival  of  fallout. 

Trade-offs  between  movement  in  fallout  and  primary  shelter  PF  are  examined 
in  a  general  manner  because  of  their  importance  in  the  transportation  algorithm 
solution. 

2.  Results 

a.  General 

Each  city  is  examined  "microscopically":  Radial  frequency  distri¬ 
butions  of  distances  from  city  centers  of  both  population  and  shelter 
are  plotted;  histograms  of  shelter  spaces  by  PF  category,  and  numbers 
and  costs  of  improvable  spaces  are  analyzed.  Improvable  space  numbers 
and  costs  average  about  17  percent  of  the  population  and  $7  per  space. 

It  was  hoped  that  the  effects  both  of  protection  factor  (PF)  and 
juxtaposition  of  people  and  shelter  and  the  latter's  effect  on  shelter 
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utilization  could  be  examined.  However,  shelter  PF's  are  comparable  in 
most  cities  (217,-317.  of  total  shelter  space  was  100  PF  or  better  for  six 
of  the  eight  cities;  the  other  two  were  407.  and  487.).  The  effects  of  dif¬ 
ferent  patterns  of  shelter  utilization  are  discussed  below. 

b .  Discussion 

For  reasons  cited  below,  the  conclusions  drawn  from  this  analysis  re¬ 
late  to  shelter  utilization  rather  than  to  the  more  fundamental  variables 
of  knowledge,  time,  and  distance  which  determine  utilization.  These 
variables  are  defined*  as  follows: 

(1)  Knowledge  by  the  population  of  what  shelter  they  should  use 
(shelter  assignment  plan,  discipline,  crowd  control,  etc.). 

(2)  Time  required  to  reach  shelter  (warning  time,  start-up  time, 
travel  time,  searching  time,  interference  of  traffic  flow,  etc.). 

(3)  Distance  to  shelter  (influences  both  time  and  knowledge,  for  the 
knowledge  will  be  less  if  distances  are  long). 

In  a  complete  analysis,  relations  among  these  three  variables  must 
be  derived  (for  obviously  they  are  not  independent),  and  these  in  turn 
must  be  related  to  the  shelter  utilization  pattern. 

In  the  analyses  reported  here,  knowledge,  time,  and  distance  to  the 
shelter  are  introduced  by  varying  "movement  restrictions"  as  follows: 

(1)  Movement  Restricted  to  Own  Standard  Location.  Knowledge  very 
great,  and  time  and  distance  allowance  great  enough,  to  allow 
best  shelter  utilization  within  one's  own  standard  location  (SL) 
(Itovement  Restriction  1). 

(2)  Movement  Restricted  to  2  Miles  from  Own  Standard  Location. 
Knowledge  very  great,  and  time  and  distance  allowance  great 
enough  to  allow  best  utilization  of  any  shelter  within  one's 
own  standard  location  (SL),  or  an  SL  not  more  than  two  miles 
away  (rectangular  distance)  (Movement  Restriction  2). 

(3)  Movement  Unrestricted  Within  SMSA.  Knowledge  very  great,  and 
time  allowance  great  enough  to  allow  best  shelter  utilization 
at  long  distances  anywhere  in  the  SMSA  (Movement  Restriction  3). 

(4)  Optimum  Shelter  Utilization  (Transportation  Algorithm).  Know¬ 
ledge  exact,  to  allow  optimal  shelter  utilization  (minimum 
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rasuaj  Lies)  anywhere  within  the  SMSA,  including  a  prior  calcula- 
l 1  an  of  dose  received  during  transit  in  fallout  and  in  the 
primary  shelter,  all  before  or  during  fallout  arrival;  speed 
is  restricted  to  2  mph .  (Movement  Restriction  A). 

An  optimum  mathematical  solution  presented  in  Movement  Restriction  4 
is  not  practically  realizable,  as  it  assumes  "perfect"  prior  knowledge  of 
the  attack  environment  and  "perfect"  population  response.  However,  a 
similar  shelter  utilization  pattern  could  be  generated  by  other  assumptions 
that  are  operationally  more  realistic;  for  example,  this  can  be  seen  for 
unrestricted  movement  (Restriction  3)  in  Table  I,  Summary  of  Casualty 
Estimates  CL  of  Population)  for  Reference  Intensity  of  2000  r/hr  and  Time 
of  Arrival  of  1,5  Hours.  Casualties  differ  between  Restrictions  3  and  4 
in  only  two  of  the  eight  SMSA's.  Unrestricted  movement  and  a  mathematically 
optimum  solution  offer  a  benchmark  of  ideality  against  which  to  measure 
other  shelter  utilization  patterns. 

TABLE  I 

Summary  of  Casualty  Estimates  CL  of  Population)  for 
Reference  Intensity  of  2000  r/hr  and  Time  of  Arrival  of  1.5  Hours 


SMSA 

Shelter/ 

Population 

Ratio 

Categories 

1  and 

Better 

Restriction  1 

Movement 

Restricted 

To  Own  SL; 

Movement 

Completed 

Before 

Fallout 

Restriction  3 

Movement 

Unrestricted 

Within  SMSA; 

Movement 

Completed 

Before 

Fallout 

Restriction  4 

Constant  Movement 
Speed  of  2  mph 
Restriction; 

Movement  Permitted 
in  Fallout  (Trans¬ 
portation  Algorithm) 

Orlando 

.20 

94 

78 

78 

Tucson 

.24 

89 

74 

74 

Charleston 

.33 

83 

70 

70 

Mob  i  1 e 

.48 

78 

54 

54 

Macon 

.50 

77 

46 

46 

Pittsfield 

.  6b 

64 

43 

43 

Waterloo 

1 . 10 

83 

0 

25 

Richmond 

1.60 

66 

0 

36 

c.  Major  Conclusions 

Efficient  shelter  utilization  appears  to  be  a  critical  problem  when 
the  shelter  spaces  exceed  population;  inefficient  utilization  can  lead  to 
a  large  number  of  avoidable  casualties.  This  is  concluded  from  the  ob¬ 
servation  of  the  rapid  increase  in  casualties  for  a  specified  attack 
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environment  as  one  restricts  movement  progressively  from  the  SMSA  down  to 
the  SL  in  which  the  shelteree  finds  himself  at  the  time  of  warning.  For 
the  most  shelter -abundant  cities  (Waterloo  ant!  Richmond  in  Table  I),  the 
percentage  of  casualties  under  a  2000  r/hr  fallout  environment  is  progres¬ 
sively  reduced.  The  percentages  reduce  progressively  from  a  two-city 
average  of  757.  to  307.  for  decreasing  movement  constraints,  and  to  07,  for 
unrestricted  movement  within  the  SMSA. 

For  the  two  most  shelter-poor  cities  (Orlando  and  Tucson),  the  cor¬ 
responding  reduction  is  from  a  two-city  average  of  917,  casualties  to  a 
two-city  average  of  767,  casualties,  a  possible  savings  of  157.  of  the  popu¬ 
lation  who  would  otherwise  be  casualties.  This  is  to  be  compared  to  the 
757.  savings  in  the  more  shelter -abundant  cities.  The  lack  of  sensitivity 
to  movement  restrictions  in  shelter-poor  cities  generally  is  due  to  the 
fact  that  there  are  enough  people  in  any  standard  location  to  completely 
utilize  available  shelter  in  that  SL.  Bius,  even  the  most  restrictive 
shelter  utilization  constraint  leads  to  good  utilization. 

Within  broad  limits,  the  higher  is  the  ratio  of  shelter  to  poDulation, 
the  greater  is  the  payoff  in  casualties  saved  by  shelter  utilization  plan¬ 
ning. 

Insofar  as  movement  distance  restrictions  are  derived  from  time 
available  for  sheltering,  one  can  deduce  the  substantial  payoff  potential 
in  effective  means  of  warning,  shelter  utilization  planning,  and  indoctrina¬ 
tion  of  the  population--particularly  in  shelter -abundant  areas. 

d.  Secondary  Conclusions 

(1)  Movement  in  fallout  to  a  better  shelter  theoretically  is  war¬ 
ranted  in  many  circumstances.  The  quantitative  trade-offs  between 
protection  factor  of  the  shelter  and  movement  time  in  fallout  are 
shown  in  Figure  31  of  Volume  H.  The  feasibility  of  applying  this 
principle  in  the  transattack  environment  is  questionable,  however. 

It  is  a  special  case  of  remedial  movement,  but  with  a  very  short 
decision  time  and  a  changing  environment. 

(2)  Neither  city  size  nor  population  density  (one  of  the  city-selec¬ 
tion  criteria),  which  might  be  supposed  to  have  an  effect  on  juxta¬ 
position  of  people  and  shelter,  could  be  shown  to  be  significant  in 
this  analysis.  This  is  concluded  from  examining  the  difference  in 
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casualties  in  each  attack  environment,  between  Movement  Restriction 
4,  (movement  in  fallout,  but  limited  movement  speed),  and  Restriction 
3  (complete  freedom)  as  a  function  of  city  size  and  population  density. 
Small  effects  which  likely  exist  either  are  swamped  by  the  more  criti¬ 
cal  shelter/population  ratio,  or  are  undistinguishable  from  the  "noise"  of 
variable  PF  distributions  in  cities  otherwise  similar  (size  and  popu¬ 
lation  density).  The  casualties  differ  between  Restrictions  3  and  4 
only  for  Waterloo  (low  density,  small  size)  and  Richmond  (high  density, 
large  size)  at  the  intermediate  attack  environments. 

e .  Recommenda  t ions 

This  analysis  shows  the  critical  sensitivity  cl  casualty  estimates 
to  shelter  utilization,  and  therefore,  to  the  corresponding  movement-to- 
shelter  assumptions,  particularly  for  shelter-abundant  areas.  Accordingly, 
further  research  is  recommended  to  give  a  more  quantitative  insight  into 
the  problem  of  movement-to-shel ter  (tine,  training  and  discipline,  and 
distance  to  move).  Better  understanding  can  lead  to  more  effective  opera¬ 
tional  plans  and  programs  in  terms  of  preattack  shelter  assignment  plan¬ 
ning,  warning  requirements,  training,  etc. 

C.  A  Generalized  Sensitivity  Analysis  of  CD  Systems 
1 .  Approach  to  the  Research 
a.  Introduction 

This  subtask  employs  specific  models  which  estimate  the  total  dose, 
Maximum  ERD,  or  probability  of  casualty  or  fatality  for  an  individual 
exposed  to  a  particular  radiological  environment.  The  total  dose  model 
is  analytical  and  the  ERD  model  is  computerized.  These  models  are  used 
to  examine  the  sensitivity  of  the  dose  (or  probability  of  casualty)  to 
variations  in  the  input  parameters  used  to  define  the  radiological  environ¬ 
ment.  The  specific  parameters  used  are  those  in  Table  II,  Input  Parameters 
for  Generalized  Sensitivity  Analysis . 

Sensitivity  indices  are  calculated  for  each  parameter.  The  sensi¬ 
tivity  index  is  defined  as  the  fractional  change  in  dose,  probability  of 
casualty  or  fatality  divided  by  the  corresponding  fractional  change  in 
the  input  variable.  Hence,  the  relative  ranking  of  these  indices  serves 
to  identify  the  input  parameters  which  influence  the  system  performance 
most  significantly. 
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b. 


ERD  Model 


In  the  computerized  ERD  model,  a  range  of  variations  significant  for 
each  input  parameter  is  estimated.  For  each  input  parameter  at  least 
three  values  are  chosen  to  cover  the  range;  these  are  listed  in  Table  II. 
Maximum  ERD  is  calculated  for  each  combination  of  input  parameters. 

These  thousands  of  data  points  then  are  used  to  obtain  an  analytical 
function  for  Maximum  ERD  which  is  equivalent  to  the  computer  model.  The 
analytical  function  allowed  the  direct  calculation  of  sensitivity  indices. 

A  similar  analysis  is  made  using  as  the  measure  of  effectiveness  the 
length  of  stay  in  primary  shelter  required  to  allow  ERD  to  fall  to  80v  or 
5  Or . 

TABLE  II 

Input  Parameters  for  Generalized  Sensitivity  Analysis 


Parameter 

Symbol 

Input  Values 

Reference  Intensity 
(r/hr) 

*o 

300,  1500,  2700,  3900 

Time  of  Arrival  of 

Fallout  (hrs) 

ta 

1,  4,  7 

Time  Outside  in  Fallout 

Before  1st.  Shelter  (hrs) 

T2 

0.0,  0.3,  0.6 

Radiation  Decay  Exponent 
(dimensionless) 

Z 

1.0,  1.2,  1.4 

ERD  Recovery  Fraction 
(dimensionless) 

F 

0.85,  0.90,  0.95 

ERD  Recovery  Rate 
(fraction/day) 

B 

0.020,  0.025,  0.030 

,  * 

Duration  of  Fallout  Buildup 
(dimensionless) 

E 

1.13,  3.64,  6.15 

Protection  Factor 
(dimensionless) 

PF 

2,  10,  40,  100,  500 

k 

E  is  a  factor  such  that  E  . 

*  duration  of 

fallout  deposition. 
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c .  Total  Dose  Model 

When  total  dose  is  used  as  the  measure  of  effectiveness  instead  of 
ERD,  it  is  convenient  and  more  economical  to  use  non-computer  analytical 
techniques.  With  a  non-computer  analytical  model,  sensitivity  indices 
are  calculated  for  all  the  parameters  listed  in  Table  II,  except  ERD 
Recovery  Fraction  (F)  and  ERD  Recovery  Rate  (B).  When  these  indices  are 
evaluated  at  the  "mean  values"  of  the  input  parameters,  they  are  found 
to  agree  closely  with  the  ERD  model.  This  non-computer  analytical  model 
is  used  to  plot  the  sensitivity  indices  as  functions  of  several  of  the 
input  variables.  These  plots  are  used  to  examine  sensitivity  indices  for 
input  parameters  at  values  other  than  "mean  values." 

Results 

a.  Introduction 

The  general  definition  of  sensitivity  stated  above  is  applied  to 
civil  defense  systems  analysis  by  relating  a  range  of  input  parameters  to 
a  range  of  uncertainty  for  physical  constants  (such  as  decay  exponent, 
bodily  repair,  as  reflected  in  the  two  ERD  parameters,  F  and  B,  etc.)  or 
to  uncertainty  in  the  attack  environment  (reference  intensity,  time  of 
arrival,  etc.).  When  this  is  done,  we  can  conclude  that  particular  care 
must  be  exercised  in  CD  systems  evaluations  relating  to  these  sensitive 
parameters,  or  that  for  parameters  to  which  probability  of  casualty  or 
time  required  in  primary  shelter  are  sensitive,  further  research  may  be 
required  to  narrow  the  range  of  uncertainty.  When  it  is  found  that  the 
output  is  insensitive,  perhaps  we  can  drop  any  variation  in  the  parameter 
in  our  shelter  system  evaluation,  or,  at  least,  we  can  be  secure  in  the 
knowledge  that  no  seriously  incorrect  results  can  arise  from  this  quarter. 

b.  Casualty  Computation  Sensitivity  Indices 

Tabl~  III,  Relative  Importance  of  Parameters  in  the  Fallout  Shelter 
Sensitivity  Analysis,  summarizes  the  primary  conclusions  of  these  sensi¬ 
tivity  analyses  which  were  concerned  with  fallout  casualty  computations. 
The  table  shows  (for  each  input  parameter  examined)  the  sensitivity  index 
obtained  using  both  the  computerized  ERD  model  and  the  total  dose  model. 
These  sensitivity  indices  indicate  the  relative  importance  of  the  parame¬ 
ters  as  they  appear  in  the  equations  used  to  compute  fallout  casualties. 
However,  this  expression  of  the  mathematical  sensitivity  of  the  parameters 
must  be  further  examined  with  respect  to  the  estimated  accuracy  or 


"uncertainty"  of  the  parameter.  For  example,  the  Radiation  Decay  Exponent 

has  a  sensitivity  index  of  -4.22,  which  indicates  that  a  +10%  change  in 

*1  2 

this  exponent  (e.g.,  a  change  in  the  radiation  decay  rate  from  t  '  to 
-1  3 

about  t  ‘  )  will  produce  a  -42.2%  change  in  Maximum  ERD  and  a  similar 
change  in  computed  fallout  casualties  for  the  appropriate  ERD  ranges. 
However,  examinations  reported  in  Volume  HI  indicate  that  errors  of  this 
magnitude  in  predicting  the  decay  exponent  are  less  likely  to  occur  than 
error  of  10%  in  predicting  other  parameters  (e.g.,  reference  intensity). 

To  account  for  such  differences  in  probable  errors  inherent  in  estimating 
the  parameter  values,  the  column  labeled  "Uncertainty"  in  Table  III  is 
employed.  The  listed  values  of  uncertainty  are  multiplied  by  the  sensi¬ 
tivity  indices  to  obtain  a  revised  ranking  of  the  parameters,  labeled 
"Approximate  Relative  Importance." 

TABLE  III 

Relative  Importance  of  Parameters  in  the  Fallout 
Shelter  Systems  Analysis 


Parameter 

Sensitivity  Index 

ERD  Model  Total  Dose 

Model 

Uncer¬ 

tainty 

Approximate 

Relative 

Importance 

Radiation  Decay  Exponent 

-4.22 

-4.36 

25% 

1.07 

Reference  Intensity 

+1.02 

+1.00 

75 

0.76 

Protection  Factor 

-0.96 

-0.99 

80 

0.78 

Time  of  Arrival  of  Fallout 

-0.50 

-0.50 

75 

0.38 

Duration  of  Fallout  Buildup 

-0.18 

-0.25 

70 

0.15 

ERD  Recovery  Rate 

-0.14 

*  - 

70 

0.10 

Time  Outside  in  Fallout 

+0.05 

+0.02 

100 

0.04 

ERD  Recovery  Fraction 

-0.30 

•  “ 

10 

0.03 

It  should  be  noted  that  the  range  of  uncertainty  is  difficult  to 
establish  for  some  of  the  parameters.  It  often  depends  upon  the  context 
of  the  systems  evaluation.  However,  Table  III  shows  that  only  minor  re¬ 
ranking  of  the  parameters  occurs  when  sensitivity  indices  are  multiplied 
by  the  range  of  uncertainty.  That  is,  the  differences  among  the  various 
"uncertainty"  estimates  are  considerably  smaller,  generally,  than  the 
differences  among  the  mathematical  sensitivity  indices. 

These  sensitivities  are  evaluated  at  the  mean  values  of  Che  parameters. 
Some  of  them  change  significantly  for  different  values  of  the  input 
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parameters.  Details  of  variations  are  contained  In  Volume  III.  The  most 
significant  arc  summarized  below. 

c.  Sensitivity  Indices  for  Time  Required  in  Primary  Shelter 

For  time  required  in  primary  shelter  to  achieve  bodily  repair  such 
that  ERD  is  50r-80r,  Che  rankings  remain  the  same.  However,  because  of 
differences  in  formulation  (no  time  outside  in  fallout)  only  reference 
intensity  (1^),  protection  factor  (PF) ,  and  radiation  decay  exponent  (Z) 
can  be  compared. 

3.  Cone lus ions 

a.  Casualty  calculations  are  quite  sensitive  to  errors  in  the  field 
decay  exponent.  This  sensitivity  remains  high  and  essentially  constant 
throughout  the  range  examined. 

b.  Sensitivity  to  variations  in  fallout  reference  intensity  and  protec¬ 
tion  factor  are  high  over  the  whole  range  of  parameter  values.  Hence,  pre¬ 
cise  knowledge  of  the  fallout  shelter  posture  and  the  fallout  reference 
intensity  is  much  more  essential  to  accurate  vulnerability  analysis  than 
the  remaining  parameters  in  most  cases  of  interest. 

c.  Sensitivity  to  time  of  arrival  of  fallout  can  be  quite  high  in  some 
radiological  environments. 

d.  Sensitivity  of  casualty  computations  to  the  remaining  parameters  is 
low  In  most  cases  of  interest. 

c.  Expansion  of  the  sensitivity  analysis  to  include  parameters  other 
than  fallout,  which  define  the  total  casualties  from  a  given  attack  on  the 
United  States,  is  necessary  before  further  conclusions  concerning  a  national 
vulnerability  analysis  can  be  drawn. 

4.  Recomncndat Ions 

a.  Vulnerability  analyses  should  employ  the  best  available  information 
on  protection  factors,  and  research  and/or  surveys  to  Improve  protection 
factor  data  should  bo  encouraged. 

b.  Because  of  the  sensitivity  of  systems  analysis  results  to  the  field 
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decay  exponent,  continuing  analysis  of  the  validity  of  the  t  *  decay 
law  should  be  made. 
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c.  Additional  atudy  la  required  to  eatabllah  the  sensitivity  of  fatalities, 
casualties,  and  dose  to  duration  of  shelter  stay. 

d.  Sensitivity  analysis  should  be  extended  to  include  the  parameters 
defining  the  effects  of  blast  and  alternative  measures  of  effectiveness 
(casualties  by  type,  dose  distribution  of  survivors,  etc.). 

e.  Sensitivity  analysis  should  be  applied  to  identifying  the  important 
cost/effectivcness  parameters  used  in  the  budget  allocation  model  (see 
Volume  I). 


